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The negatively charged nitrogen-vacancy (NV) center spin in diamond can be used to realize 
quantum computation and to sense magnetic fields. Its spin triplet consists of three levels labeled 
with its spin z-components of -|-1, 0, and -1. Without external field, the -fl and -1 states are 
degenerate and higher than the 0 state due to the zero-field splitting. By taking the symmetrical 
and anti-symmetrical superpositions of the -fl and -1 states as our qubit basis, we obtain exact 
evolution operator of the NV center spin under time-dependent magnetic field by mapping the three- 
level system on time-dependent quantum two-level systems with exact analytical solutions. With 
our exact evolution operator of the NV center spin including three levels, we show that arbitrary 
qubits can be prepared from the starting 0 state and arbitrary rapid quantum logic gates of these 
qubits can be realized with magnetic fields. In addition, it is made clear that the typical quantum 
logic gates can be accomplished within a few nanoseconds and the fidelity can be very high because 
only magnetic field strength needs to be controlled in this approach. These results should be useful 
to realizing quantum computing with the NV center spin systems in diamond and exploring other 
effects and applications. 

PACS numbers: 75.75.-c, 03.67.-a, 75.10.-b, 75.90.+W 


I. INTRODUCTION 

The negatively charged nitrogen-vacancy (NV) cen¬ 
ter in diamond has been intensively investigated because 
it can be used for realizing quantum computation and 
sensing weak magnetic field, electric field, strain etcii^— 
As quantum technology evolves, one can manipulate the 
NV center spin with electromagnetic field^^— , optical 
fielder—, and stress fiel d^^d^ . Its spin ground state is 
a triplet {S = 1) with 0 and ±1 as the spin z exponents, 
and its zero-field splitting makes the |0) state the lowest, 
leaving the | + 1) and | — 1) degenerate. One can apply 
a longitudinal magnetic field to split the | + 1) and | — 1) 
states, adding a Zeeman term to the spin Hamiltonian. 
Usually, it is preferable to work with a two level system 
representing a qubit by using the |0) state and one of the 
I + 1) and I — 1) states, and then control it with microwave 
pulse or radio-frequency wave^. In such approaches, how¬ 
ever, one usually use only weak driving under the rotat¬ 
ing wave approximation, or else the third state cannot 
be neglected. Additionally, one needs accurate timing to 
track the relative phase between the |0) and | -I- 1) (or 
I -I- 1)) states due to the zero-field splitting effect. Re¬ 
cently, the NV center spin was used to experimentally 
realize geometric quantum gates with high fidelit y^^d^ . 
In this approach, the | + 1) and | — 1) states are used to 
encode the qubit. The latest report shows that the | + 1) 
and I — 1) states can be manipulated with internal strain 
fieldi^. Nevertheless, the geometric quantum gates need 
half microsecond to be accomplished and the mechani¬ 
cal quantum gates require longer time. Therefore, much 
more rapid quantum logic gates of such qubits are highly 
desirable. 

In this paper, we use |±) = (| -|- 1) ± | — l))/-\/2 of the 


NV center spin as the qubit basis and then show how 
to prepare arbitrary qubits from the starting |0) state 
and how to construct rapid quantum logic gates for the 
qubits. On the basis of powerful exact analytical results 
of a quantum two-level system under a time-dependent 
magnetic fiel d^^i^^ , we construct exact evolution opera¬ 
tor for the NV center spin (including three levels) under 
time-dependent magnetic fields. Then, we use these ex¬ 
act evolution operators to prepare arbitrary qubit states 
with the basis |±) from the starting state |0) and realize 
arbitrary quantum logic gates. All the typical quantum 
logic gates can be completed within a few nanoseconds. 
The fidelity can be made very high because one needs 
to control magnetic field strength only. More detailed 
results will be presented in the following. 

The rest of the paper is organized as follows. In Set. 
II, we define the Hamiltonian and elucidate the new spin 
basis. In Sec. HI, we construct exact evolution operators 
for the NV center spin under time-dependent magnetic 
fields. In Sec. IV, we show how to use the exact evolu¬ 
tion operators to prepare arbitrary qubit states with the 
basis |±) from the starting state |0). In Sec. V, we con¬ 
struct arbitrary quantum logic gates for the NV center 
spin qubits by using the exact evolution operators. Fi¬ 
nally, we make some necessary discussions and give our 
conclusion in Sec. VI. 

II. NEW SPIN BASIS 

The Hamiltonian of NV center spin S, in 
the presence of time dependent magnetic field 
B{t) = {Bx{t),By{t),Bz{t)), can be written as 

H = +jS ■ B(t), 


( 1 ) 


2 


where h = 1 is used, D = 2.87GHz is the zero-field split¬ 
ting, and 7 = 2.8MHz/G is the electron gyromagnetic 
ratio. Accordingly, the Schrddinger equation for time- 
evolution operator U is given by i-^U = HU. 

In the Sz representation, the matrix form of the Hamil¬ 
tonian o can be written as: 


where J{t)=jBi{t). Because Ht consists of a 2x2 block 

= (m T) («) 

and a 1x1 block Hti = D, we can focus on Ht 2 in Eq. 
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Introducing a 2x2 unitary transform 
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1 
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Introducing the unitary transform 



= 0 1 0 
\ J_ 0_^ 

\ V2 V2 
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we can change Ht 2 into Hh 2 = H 2 + 7?/2, with H 2 given 

by 


H 2 




( 10 ) 


we can transform the Hamiltonian (I2|) into 


D 

iBx 

iBz \ 


iBx 

0 
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1 

to 
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It means that the basis is changed from (| -|- 1),|0),| —1)) 
to (|+),|0),|-)), where |+) = ; 75 (| -k 1) -k | - 1)) and 
|—)=-^(|-kl) — I — 1)). When the magnetic field is turned 
off, the D term will produce the same overall phase for 
both |-k) and |—). Therefore, there will be no relative 
phase between |-k) and |—), and |-k) and |—) can be used 
to make a stable qubit basis. 


The constant term in Hh 2 contributes only an over¬ 
all phase to the time evolution operator. For the time- 
dependent two-level Hamiltonian H 2 , there are many 
exact solutions such as thoseii^— . Here, we choose a 
powerful method to construct its exact time evolution 

operator-i2i2i. 

The Schrddinger equation of H 2 can be exactly solved 
by the evolution operato r^^^^ , 

4 r) 

and the matrix elements un and U 21 and the quantity 
J{t) can be expressed 


III. EXACT EVOLUTION OPERATORS 

With a special magnetic field B{t)={aBi, PBi,0), the 
Hamiltonian in Eq. o becomes 

/ D ayBi 0 \ 

Ho = ayBi 0 ijSyBi , (5) 

\ 0 -iPyBi D J 

where a and j3 are adjustable real parameters satisfying 
o? -k /3^=1. Because of this condition, a and (3 can be 
parameterized as a = cos0 and /3 = sin0, where —tt < 
6 < TV. Introducing another unitary transform 

( COS0 0 isin0 \ 

010 (6) 
isin0 0 COS0 J 

we can transform Ho into a block-diagonal matrix Ht = 

S2H0SI 


uii{t) = cos(x(t))e*«-b) 

U2i{t) = ir7sin(x(t))e®^+^‘) 

^W = 7rf|^-5v/^D^^4Wcot(2x(t)) ’ 

( 12 ) 

where is defined as 

i± = j ^ - 4x2 CSC (2x) ± ^ arcsin( ^ ) ± 7 ^ . 

(13) 

Here, r] can take either -kl or —1, and x must satisfies 
three conditions: |x(t)| < -j, x(0) = 0, and X(0) = —? 7 -j. 
By choosing suitable x(0j can exactly construct J{t) 
in H 2 and the evolution operator U 2 in this way. 

After adding the phase factor due to the T>/2 term and 
making the inverse unitary transformation with T^, we 
obtain the evolution operator for Ht 2 , 


where the matrix elements un and U 21 are expressed as 


/ ^ 

j{t) 

0 \ 




J(t) 

0 


(7) 

f Unit) 

= cos(x(t))cos(^_(t)) -ki??cos(^+(t))sin(x(t)) 

V 0 

0 



\ U2l{t) 

= ? 7 sin(x(t))sin(^+(t)) -kicos(x(t))sin(^_(t)) 
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Consequently, the whole time evolution operator of the 
Hamiltonian Hf can be written as: 

/ Mil -U21 0 \ 

Ut (t) = U21 m!i 0 (15) 

\ 0 0 / 

After making the inverse transform we can get the 
time evolution operator of the starting Hamiltonian Ho 
in the new basis of |+) and |—): 

Uoid, t) = d{t)x 

uiio^ + d{t)j3‘^ —U2ia —i(d(t) — uii)aid 

■U 21 Q: u\i *M2i/? 

i [d(t) — uii) a(3 iu 2 iP d(t)a'^ + un/d'^ 

where a = cos0, /3 = sin0, and d{t) = 

From above equations (fT^ and (IT^ , we can see that if 
the time-dependent function x(t) is specified, the evolu¬ 
tion operator (1161) will be determined immediately. Using 
this time evolution operators, we can control a single NV 
center spin exactly and efficiently, and thus we can ini¬ 
tialize and gate arbitrary qubits in the stable basis of |-|-) 
and |—). 



IV. EXACT INITIALIZING OF ARBITRARY 
QUBITS 


Experimentally, the NV center spin can be easily pre¬ 
pared in state |0). We try to realize state transfer be¬ 
tween |0) and |±). With the time evolution operator Uo 
applied, the state |0) will become 


Uo{e,t)\Q) = d{t) 


— cos Oirj sin x sin C-i- ~ * cos x sin ) 
cos X cos — irj sin x cos 
i sin Oiji sin x sin 'C+ — i cos x sin ) 

(17) 


We require that the function x(i) is given by 


X it) = At 


2KX^t^ 

3 




(18) 


where A is defined as -jj k is an adjustable parameter, 
and Tf describes the time duration. Using r]=—l, we 
have ^+(T/)=^_(T/) according to equation (IT^ . Because 
the target state doesn’t contain state |0), we need to set 
uii{Tf) = 0, i.e. cos(^+(T/)) = cos(^_(r/)) = 0. Then 
the quantity xi^f) contributes an overall phase in the 
state UoiO,Tf)\Q) in Eq. (fTTl) . In order to achieve a 
minimal time value Tj and a finite field pulse in the time 
interval t € (0, Tf), we need two conditions: 0 < x(^/) ^ 
^ and 


J dt^y/D'^ - 4x2 csc(2x) = (19) 


Once we set t = Tf and choose a value for xi^f), 
can be formally solved by using equation (1181) . reading 


K = 15{XTf — xiTf))/(XTf)^. The time duration Tf can 
be solved from equation dT51) . and then k can be calcu¬ 
lated immediately. In the following, we shall show how 
to initialize three typical qubits from the spin state |0). 

Initializing the basis states |±). Choosing a = 0 (or 
/3 = 0) in Eq. (5), we can easily get the target state 
C7o(|,r/)|0) = |-) (or Uo{0,Tf)\Q) = |-|-)) with an over¬ 
all phase. In this way, we get |±) from |0). 

Initializing a superposed state cos6(i|-|-) -|- isin0i|—). 
In this cases, we can assume 0 < < tt without losing 

any effective information. Choosing a and (3 in Eq. (5) 
to satisfy the equality arctan(/3/Q;) = 0 > 0, we can let 
9 = TT — 9i in Eq. (17) and thus obtain the final state 

( cos 01 \ 

0 (20) 
i sin 01 j 

Neglecting the overall factor, we obtain cos0i|-|-) -I- 
*sin0i|—). 

Initializing arbitrary state cos 0i|-|-)-I-sin 0i|—). We 
assume that 0 < p <Tr. State cos0i|-|-) -|-tsin0i|—) can 
be obtained from |0) in the same way as above. The next 
step is to realize the phase factor for |—). If Bi{t) 

in Eq. ([5]) is a time-independent magnetic field Bq, the 
method presented in Sec. HI is not applicable, but it is 
not difficult to derive the evolution operator: 

P{t, ao,/3o) = d(T)x 

Podir) + (-l)"ao 0 iaoPo ((-!)" - d(r)) 

0 (-1)" 0 
-iaoPo ((-!)" - d('r)) 0 ald{T) + (-l)"/3o 

(21) 

where n is a non-negative integer. We assume T=Tf when 
the evolution ends. Here, t/ must satisfy the equality: 

-jT/Y^T-i-4^^2^=ii7r. If n = 0, we must have Tf = 0. 
Consequently, it means g) — ^ = 0, and the operator 
P(t/, ao)/3o) becomes identity operator, or trivial op¬ 
erator. For any nontrivial evolution operator, we have 
Tf > 0 and n = 1, Because of ^ G ^], we need 
nontrivial evolution operators for |<p — >0. It is easily 

seen that only if ao/^o = 0 in (21), can arbitrary relative 
phases be realized by the operator P(r/, ao,/3o)- This 
means that the magnetic field should be applied in either 
x-axis or y-axis. For ip £ [0, !■), we need to set ao = 0 and 

then obtain Tf = ^{(p + ^) a.nd Bq = - 1; 

and for G (^, tt], we need to set /3o = 0 and then obtain 
= %i^ - t) and Po = ^ We can pa¬ 

rameterize ao = cos 00 and /3o = sin 0o. In this way, time 
duration r/ will be only few nanoseconds and magnetic 
field Bq will be reasonable at the same time. Because 0o 
(ofo and /3o), Tf, and Bq are determined by p, we can use 
Pf{(p — to denote the evolution operator. 

Therefore, the whole procedure, achieved in two steps, 
can be represented as the evolution operator 

Ui{g,,9,)=Pf{g>-^)Uoi7r-9,,Tf) (22) 
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FIG. 1. The final magnetic field B\{Tj) in unit of 10®G de¬ 
pending on parameter x{Tf)- 


For practical application, it is useful to adjust the value 
x{Tf) to connect the two magnetic fields Bi{t) and Bq 
at the time Tf along either x-axis or y-axis. We show in 
Fig. 1 that magnetic field Bi{Tf) can vary from 450G 
to 740000G when x{Tf) changes within (0,^]. Because 
of the large domain of Bi{Tf), we can likely realize con¬ 
tinuous connection of magnetic field in either x-axis or 
y-axis. 


V. REALIZATION OF QUANTUM LOGIC 
GATES 

In last section, we show how to initiate the basis states 
|±) and their superposed states from the state |0). Our 
arbitrary qubits are made from the basis states |±). In 
the quantum circuit model of computation, a quantum 
gate is a basic quantum circuit operating on a small num¬ 
ber of qubits. We show how to realize typical quantum 
gates on the qubit in the following. In some of the cases, 
the special state |0) can be used as a auxiliary state, or 
a bridge. 

j phase shift gate: |-|-) -|- |—) —|-|-) -|- |—). In this 
case, we need only a phase factor e^^ for the |—) term. 
This can be achieved by applying a unitary transforma¬ 
tion F’/(f) on the starting state |-|-) -|- |—). As a result, 
we obtain the evolution operator for the j phase shift 
gate: 

U^=Pf{^)e-^\ (23) 

The time duration ri and the time-independent magnetic 
field hi can be given by ti = ^ ^ 1.1ns and hi = ^ 

888G. ^ 

I phase shift gate: |-|-) -|- |—) —>• |-|-) -|- e*T|—). It is 
similar to the j phase shift gate. The phase factor is 
in this case. Applying Pf{j) on (|-|-) -|- |—)), we obtain 
the evolution operator for the ^ phase shift gate 

C/. =P/(J)e-*U (24) 


The time duration T 2 and the time-independent magnetic 
field h 2 are given by T 2 = ~ 1.64ns and /12 = ^ 

452G. 

Pauli-X gate: |-|-) —>■ |—). This gate can be realized 
by applying {7j(0,ti) on the initial state |-|-) and then 
Go(f ,^ 2 ) on the resulting intermediate state | 0 ): 

|-)cxC/„(|,t2)C/t(0,G)|+). (25) 

In this way the gate is realized in two steps. Letting 
Xi(^i) = X 2 (^ 2 ) = f, we have G = ^2 = T, and then 
obtain ki = K 2 = Using the phase condition 

in (fTOl) . we can obtain T = 2^. For the first step, the 
time-dependent field function J{t)/X and the probability 
^ 1 +) of |-|-) as functions of time (At) are shown in Fig. 2. 
For the second step, we have similar time dependence for 
the field and the probability. After these two steps, the 
evolution operator for the Pauli-X gate is given by 

Vx = U,i^,T)Ul{ 0 ,T)P^/^. (26) 

And the total time interval is equivalent to 2T ~ 2.7ns. 

In addition, this state can be realized without apply¬ 
ing Uo{ 02 ,Ti) on state | 0 ), where Ti is pulse time dura¬ 
tion. At this time, we must guarantee matrix element 
U 2 i{Ti)= 0 . After applying operator Uo{02,Ti) on state 
|-|-), we derive 

/ sin^ 6*2 — eos^ 02 

Uoi02,Ti)\+) = d{Ti) I 0 

\ -|-j eos 02 sin02 

(27) 

with ^±(Tj)=7r + 2rmr and m is positive integer. The 
condition to achieve state |—) is sin^( 202 )[l -|-cos(-jri — 
X3 (Ti))]= 2. It can be satisfied by setting cos(-jTi — 
X3(7’i))=l and 02=f. A reasonable result is -jTi — 
X3(7’i)=27r and then ^ 3 =^^. Using the phase condi¬ 
tion about 'f±(T'i) in equation (ITOl) by replacing | with 
TT -|- 2m7r, we can get a self-consistent value, Ti ~ 5.3ns 
from ~ 7.542 with m = 4. Therefore, the evolution 
operator can be expressed as 

[/^ = C/o(02,ri)e*(-^^^+’"/^ (28) 

The smooth pulse and probability evolution of state |-|-) 
is shown in Fig. 3. Comparing it with Fig. 2, we can see 
that it is not as efficient as that using the intermediate 
state | 0 ) because the magnetic field as a function of time 
is irregular. 

Hadamard gate: |±) —:^l+) In this situa¬ 

tion, we can’t use operator (1^^ directly. So in first step, 
we change it to state |0). It is the same step mentioned in 
dealing with Pauli-X gate: f7j(0, T)|-|-). After this step, 
for case | 0 ) turning to -^( 1 +) + |—)), we only have to 
set 61 = j and (^ = 0 in operator (1^ . The evolution 
operator for the Hadamard gate is 

G^ = c/z(0,J)c/t(0,r). 


(29) 
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FIG. 2. The time t (in 1/A) dependence of J{t)/\ (upper 

- in Eq. 


panel) and ^l+> (lower panel), with n = 
(18). 


f 



0 1 2 3 4 5 6 7 

At 

FIG. 3. The time t (in 1/A) dependence of J{t)/\ (upper 
panel) and ^1+) (lower panel), with At = in Eq. (18). 


The total time is 2T + t/ ~ ^ 3.83ns. 

In addition, this gate can be realized without using 
|0). In this case, the static field is applied in z-axis, and 
the time evolution operator for the Hamiltonian (|T]) with 
magnetic field ( 0 , 0 , B^) is given by 

( cos{Jzt) 0 —isin(J 2 t) \ 

0 6 *-°* 0 (30) 

—isiii{Jzt) 0 cos( Jzf) J 


where Jz=jBz. If the initial state is |+) and Jztf=j, 
where t/ is time duration, Q{tf) makes the state |+) 
become (5(t/)I+) = -^(1+) —*|—))■ In order to get target 

state -^( 1 +) + |—)), we need a phase factor for the 
|—). It can be shown that the time value and magnetic 
field is the same as ri and hi. With the phase shift 
operator the whole process can be represented as 


u'h = Pf{^)Q{tf) 



0 

1 


V2 

0 

i^iDtf 

0 

jj_ 

0 

1 


V2 


o-iDtf 


(31) 


The total time duration is t/ + ti = ~ 1.375ns 

with Bz ^ lOOOG. Furthermore, this gate can change |—) 
into ^|+)-^|-). 

Arbitrary gating. In this case, we need to realize 
cos 6 »i|+) + sin 6 li|-) —^ cos 6 » 2 |+) + 6 ®“^= sinezh). 
Assuming 0 < 0 i ,02 < and 0 < y:>i,^i < tt. If 6 i=02, 
we need only to modify the relative phase. According 
to the previous subsection, we can realize this control 
through the intermediate state |0). The first part is actu¬ 
ally the inverse process of that in Sec. IV. As the operator 
( 1 ^ is unitary, the first process can be achieved by the 
operator Ul{ipi,0i), with T 2 being the pulse time dura¬ 
tion. Then, the target state can be realized by Ui{(fi2,02) 
with T 3 being the pulse time duration. The whole process 
can be represented by 


Ua{‘P2,02,‘P1,0i) — Ul{(p2,02)uj{(pi,0i) (32) 


This gate can be realized without the intermediate 
state |0). It is clear from operator (1501) that Jzt can 
be seen as one variable because the overall phase fac¬ 
tor does not play role in the gate. As a result, the 
larger the longitudinal field, the smaller the time dura¬ 
tion. The initial relative phase (pi should be adjusted 
before the probability amplitude is changed. At first, 
we change the relative phase ipi into | by Pj{(pi — f), 
making a state |d’(0)) = cos0i|-|-) -I- isin0i|—). Then, 
we apply the time evolution operator Q{te) on |$( 0 )), 
where tg is the time duration, and obtain the state 
|$(te)) = (cos( 6 »i - Jzte)\+) +ism{0i - Jzte)\-))e~^^^‘. 
Here, can be expressed as 


r 01-02 if 01 > 02 

{ 01 - 02 + 271 if 01 < 02 


(33) 


Finally, we apply the phase regulation Pf{ip 2 — f) to get 
the phase (p 2 . The whole procedure can be represented 
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as the unitary operator: 

U'A{ip2,d2,<Pi,0i) = Pf{^2 - ^)Q{te)P}{<Pi - |) (34) 


TABLE I. The time durations (To) of the \ phase, ^ phase, 
Pauli-X, and Hadamard gates 


Gate 

J phase 

^ phase 

Pauli-X 

Hadamard 

Tg 

1.1 ns 

1.6 ns 

2.7 ns 

3.8 ns 




5.3 ns 

1.4 ns 


VI. DISCUSSION AND CONCLUSION 

The time durations Tq of the typical gates are sum¬ 
marized in Table 1. It is clear that the gates need at most 
a few nanoseconds. Furthermore, it can be estimated to 
take only nanoseconds to complete the initialization and 
gating of arbitrary qubits. Therefore, the qubits on the 
basis of |-|-) and |—) can be fast initialized and gated with 
magnetic fields, in comparison to those in terms of | -h 1) 
and I — 1) with strong mechanical drivingi^. When the 
magnetic field is applied along the z axis, the effect of 
the time evolution operator (1501) on our NV center qubit 
with the basis of |-|-) and |—) looks like that on spin-i 
qubit in “bang-bang” approach^^, but they are different 
from each other. Because here one needs to control the 


magnetic held strength only, the theoretical hdelity can 
be very high for these gates. 

In summary, by choosing |±), dehned as (| -I- 1) ± | — 
l))/-\/2, of the NV center spin as the qubit basis, we 
have obtained exact evolution operator of the NV cen¬ 
ter spin under time-dependent magnetic Held by map¬ 
ping the three-level system of the NV center spin on a 
two-level spin system under a time-dependent magnetic 
Held and using the existing exact analytical results of the 
quantum two-level systemi2i2i. With the exact evolution 
operator of the NV center spin including three levels, we 
have shown how to pre pare arbitrary qubits with the ba¬ 
sis |±) from the starting |0) state and to realize arbitrary 
rapid quantum logic gates for the qubits. It also has been 
estimated that the typical quantum logic gates can be ac¬ 
complished within a few nanoseconds. We believe that 
the Hdelity can be very high because only magnetic Held 
strength needs to be controlled, these results are useful 
to realizing quantum computing with the NV center spin 
systems in diamond and exploring other quantum effects 
and applications. 
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